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Abstract 

We consider spatially two dimensional Madelung fluid whose irrotational motion reduces into 
the Schrodinger equation for a single free particle. In this respect, we regard the former as a direct 
generalization of the latter, allowing a rotational quantum flow. We then ask for the most likely 
wave function possessing a given average energy by maximizing the Shannon information entropy 
over the quantum probability density. We show that there exists a class of solutions in which the 
wave function is self-trapped, rotationally symmetric, spatially localized with finite support, and 
spinning around its center, yet stationary. The stationarity comes from the balance between the 
attractive quantum force field of a trapping quantum potential generated by quantum probability 
density and the repulsive centrifugal force of a rotating velocity vector field. We further show 
that there is a limiting case where the wave function is non-spinning and yet still stationary. This 
special state turns out to be the lowest stationary state of the ordinary Schrodinger equation for a 
particle in a cylindrical tube classical potential. 

PACS numbers: 03.65. Ge; 03.65.Ta; 05.65.+b 
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I. MADELUNG FLUID: GENERALIZED SCHRODINGER EQUATION 



Let us consider a quantum system of a single free particle with mass m whose dynamics 
is confined in two dimensional space, q = {x, y}. The complex- valued wave function ^(q; t), 
where t is time, is governed by the linear Schrodinger equation 

»^(q;t) = -^^(q;t)> W 

where <9^ = <9 q • <9 q is two dimensional Laplace operator. 

Let us project the above dynamics onto real-space. To do this, let us write the wave 
function in polar form, ip(q;t) = -R(q; t) exp(zS'(q; t)/h), where the quantum amplitude 
R and the quantum phase S are real-valued functions. Inserting this into Eq. ([1]) and 
separating into the real and imaginary parts, one gets 

8 t S+ a * S ' d * S +U = 0, 
2m 

dtp + d q • = o, (2) 

where U is the so-called quantum potential generated by the quantum probability density 

p = ipip* = R 2 as 

Next, let us postulate a velocity flow generated by the quantum phase as follows 

v(q;t) = -9 q S'(q;t). (4) 
m 

Using this, the coupled equations in <^ can be rewritten into the following form 

m Tt = - d « u > 

d tP + 9 q ■ (vp) = 0. (5) 

The upper equation can then be seen as an Euler equation for the velocity flow v(q; t) 
dragged by a special kind of force field generated by the quantum amplitude, —dqU, later on 
to be referred to as quantum force. On the other hand, the lower equation can be regarded 
as the continuity equation for the quantum probability density, p(q;i). We have thus a 
non-linear fluid dynamics in real-space. This real-space hydrodynamical interpretation of 
Schrodinger equation is introduced by Madelung as soon as 1926 Moreover, the same 
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mathematical formalism is also used to based what is later called as pilot-wave interpretation, 



IB, 3- 



initiated by de Broglie and expanded by Bohm 

If we assume that the quantum phase S is twice differentiate, then the velocity flow 
defined in Eq. (J3J) is irrotational 

l/d 2 S <9 2 S\ n 
<9 q xv = - — --- - z = 0, (6) 



m\d x d y d y d x , 

where z is a unit vector orthogonal to the xy— plane. In this present work, we shall follow 
the suggestion made by Takabayashi [6] to consider the pair of equations in (jSJ) as a direct 
generalization of Schrodinger equation by allowing rotational flow. The special case of 
Madelung fluid with only irrotational flow will then reduce into the Schrodinger equation of 
(CEJ). In other words, we shall consider Eqs. (jSJ) as describing the quantum system of a single 
free particle, allowing a class of wave functions whose quantum phase possesses singularity 
such that Eq. df)J) is no more valid. 

One of the important feature of the above Madelung fluid dynamics is that it is self- 
referential. One can see that the quantum probability density p(q; t) will generate quantum 
potential U (q; t) through Eq. (jHJ), and in turn U (q; t) will tell p(q; t) the way it must evolve 
through Eqs. (j5J), and so on and so forth. There is a dynamics circularity between the 
object to be ruled, p, and the rule, U. Self-referential property is a sympton of complex non- 
linear systems, and is argued to be the general origin for the emergence of many interesting 
phenomena observed in Nature [7|. Hence, it is of great interest to evaluate the fixed point 
of the self-referential Madelung fluid dynamics. 

In this paper, we would like to address the following simple yet fundamental question: 
What is the most likely wave function of a single free particle with the average quantum 
mechanical energy (E) given by 

(E) = J ^*(q)(-|^)^(q), (7) 

where dq = dxdy. We shall show that there exist a class of solutions which turns out to be 
the fixed point of the self-referential Madelung fluid dynamics described above. 

II. A CLASS OF SELF-TRAPPED QUANTUM PROBABILITY DENSITIES 



To investigate the question set up at the end of the previous section, let us discuss 
the following readily executable problem. First, in quantum mechanics p(q) provides the 
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FIG. 1: The profile of a self-trapped quantum probability density (upper) and the corresponding 
quantum potential it generates (lower). 



essential information about the position of the particle 



101 ] . Let us quantify this 



information using differential entropy or Shannon information entropy for continuous random 



variable as 
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H[ P ] 



dqp(q) lnp(q). 



(8) 



It measures the degree of spatial localization of the quantum probability density, p(q). Now, 
let us search for a class of wave functions which maximizes the above Shannon information 
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entropy, given a finite value of the average quantum potential 

U= [ dqU(q)p(q). (9) 

way 



This is the so-called maximum entropy principle [13j . It has been argued as the only 



to infer from an incomplete information which does not lead to logical inconsistency 
Hence, it will give us the most likely p(q) whose average quantum potential is equal to U. 
Maximizing Eq. ([8]) with constraint of Eq. ([9]), one directly obtains the following relation 

p(q;t) = -^exp (-/?[%;*)), (10) 

where (3 is a real- valued constant (Lagrange constant), and Z is a normalization factor given 
by 

Z = [ dqexp(-f3U). (11) 



Notice that Eq. ( flOl) together with the definition of quantum potential given in Eq. ([3]) 
comprise a differential equation for p(q) or U(q), subjected to the condition that p(q) must 
be normalized. (3 can then be calculated by inserting Eq. ffTUj) back into Eq. ([9]). Hence, 
it is a function of the average quantum potential, (3 = (3{U). Below, instead of using U, we 
shall consider (3 as parameter of the class of quantum probability density p(q; f3) maximizing 
Shannon information entropy. Furthermore, for later convenient, we shall limit ourselves to 
positive definite f3. 

Let us show that there are infinitely many quantum probability densities satisfying Eq. 
(|T0|) . To do this, applying Laplace operator to both sides of Eq. (jTOl) . using the fact that 
p = R 2 and the definition of quantum potential given in Eq. ([3]) one obtains the following 
non-linear partial differential equation 

%U = ^U-W + ^U. (12) 

Let us remark first that the above differential equation is invariant under the following 
rotation of coordinate system 

x' = x cos 9 — y sin 9, 

y' = x sin9 + y cos9, (13) 
where 9 is the angle of rotation. Hence, if U(q) is a solution then so is U(q'). 
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We shall resort to numerical methods to solve Eq. ({12]) . For simplicity, let us confine 
ourselves to a class of solutions in which the quantum probability density is separable 

p(q) =Px(x)p y (y), (14) 

where pi(i) depends only on i = (x, y). In this case, the quantum amplitude is also separable 
i?(q) = R x (x)Ry(y) such that the quantum potential is decomposable as 

U(q) = U x (x) + U y (y), (15) 

where each term on the right hand side is given by 

= -stst- t = x ' y - (16) 

Using Eq. (fT5|) . Eq. (fl2l) can be collected as 

dlU X - \{d X XJ X f - A 2 U X = -diUy + ^(dyUyf + A^Uy. ( 1 7) 

Hence, both sides must be equal to a constant, say E. Let us choose the case when E = 0. 
One therefore obtains the following two de-coupled ordinary differential equations 

d 2 x U x = ^(d x U x ) 2 + A 2 U X , 

9 2 yUy = ^(dyUyf + A'Uy. (18) 

Fig. [T] shows the numerical solution of Eqs. fllSp with the boundary conditions d x ll x (0) = 
d y Uy(0) = and U x (0) = U y (0) = 1, for (3 = 1. For convenient, all numerical calculations 
in this paper is done by putting m = h = 1. The profile of quantum probability density 
(upper) is plotted together with the profile of quantum potential it itself generates (lower). 
We can see clearly from the numerical solution that globally the quantum probability density 
is being trapped by its own quantum potential. 

The global self-trapping property can be justified for any positive value of £/j(0) > as 
follows. First, from the upper part of Eq. fflBl . applying the boundary conditions at x = 0, 
one gets d x U x {$) = AU x (0) > 0, such that U x (x) is locally convex at x = 0. Hence, since 
d x U x (0) = 0, at spatial points nearby x = one has U x (x) > U x (0) > 0. Moreover, since the 
first term on the right hand side is always non- negative, at this region one has d 2 U x > 0. 
This geometrical reasoning can be extended such that for the whole spatial points one gets 

d 2 x U x > 0. (19) 



The same thing applies for U y . On the other hand, from Eq. (fl5|) one has d x d y U = 
0. In general one therefore obtains didjU > 0. Hence, the case when Ui(0) > will 
give an everywhere positive and convex quantum potential, q = turns out to be the 
global minimum of U(q). One can thus conclude that quantum probability density, p(q), is 
being trapped by its own self-generated quantum potential, U(q). Moreover, the rotational 
invariant of the differential equation (Tl2l guarantees that one can generate a new solution 
by rotating the solution given in Fig. [TJ We have thus a class of infinitely many self-trapped 
quantum probability densities characterized by (3. 

Next, using the fact that Ui(i) is positive everywhere and Eq. ( fl6l) . one gets 

<9^<0, i = (x,y). (20) 

Hence, Ri{i) is everywhere concave. Since Ri(i) is finite and possesses symmetrical property 
Ri(i) = Ri(—i), then Ri(i) must cross the i— axis at finite value of i = ±z m . In other words, 
Ri{i) possesses only finite support on i— axis, namely [— i m ,i m ]- One can thus conclude 
that -R(q) has only a finite support on xy— plane, that is all points q belonging to the 
rectangle [—x m , x m ]® [— y m , y m }. At the boundary of this rectangle, namely the lines i = ±z m , 
i = (x, y), U(x, y) is infinite, yet R(x, y) is vanishing. We shall show later that for a subclass 
of solutions, U remains finite for non- vanishing (3. Moreover, at % = ±i m , diRi is not 
continuous but diPi = IRidiRi is continuous and equal to zero. 

We have thus developed a class of quantum probability densities satisfying Eq. ffTO]) 
which is self-trapped by its own self-generated quantum potential. Yet, we have not specified 
the quantum phase or the velocity vector field, to completely identify our Madelung fluid 
dynamics. Nevertheless, one can see that due to the global trapping property of the quantum 
potential, any initial velocity vector field will be dragged by the quantum force to localize 
even further. In the next section, we shall specify a specific velocity vector field which 
generates a force that exactly cancels the attractive quantum force of the trapping quantum 
potential to create a stationary Madelung flow. 

III. A CLASS OF STATIONARY-SPINNING WAVE FUNCTIONS 

Let us find a class of self-trapped wave functions satisfying the differential equation ([12]) 
which is rotationally symmetric. They are the eigenfunctions of the rotation operator. 
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FIG. 2: The radial profile of a self-trapped quantum probability density (solid line) trapped by 
its own quantum potential (dashed line). The quantum potential is shifted down such that its 
minimum is equal to zero and the quantum probability density is rescaled by a constant factor. 
See text for detail. 

To do this, it is convenient to use polar coordinate {r, 6}, where r = a/x 2 + y 2 and 6 = 
tan _1 (?//x). For a rotationally symmetric solution, one imposes dgU = 0. This condition 
leads to {d x ,d y } = {cos9d r ,sm9d r }. Inserting these into Eq. ( fl2l) . one finally has to solve 

d?U+ 2 -d r U-^d r Uf-^U = 0. (21) 

Fig. [2] shows the numerical solutions of Eq. (j2*T|) for (3=1. p(r) (solid line) is plotted 
together with the corresponding U(r) (dashed line). We can see clearly that globally the 
quantum probability density is indeed being trapped by the quantum potential it itself gen- 
erates. Moreover, for a rotationally symmetric quantum potential, the attractive quantum 
force vector field, —d^U, is directing toward the origin, q = 0. 

Numerical solutions of the differential equation (j2"Tj) show that again the quantum prob- 
ability density possesses a finite support, M.. Since the quantum probability is rotationally 
symmetric, then the support is a disk with finite radius. Let us denotes the radius as r = r m . 
At the boundary of the support, dA4, which makes a circle x 2 + y 2 = r 2 ^, the quantum po- 
tential is infinite, U{r m ) = oo. See Fig. [2j In Fig. [3j we plot the radius of the support 
r m (P) and the second moment defined as r 2 (j3) = J dqr 2 p{r) against the value of (3. Both 
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FIG. 3: r m (upper) and r 2 (lower) against (3. 

are evaluated numerically. We can see that both quantities behave almost similarly as the 
function of (3. First, both increase very quickly for small (3 and then both seem to converge 
toward certain finite values for infinite value of f3 

lim r m ee r^, lim r 2 = r^. (22) 

In the next section, we shall discuss in more detail the behavior of p(r; 0) and U(r;/3) for 
oo. 

For the opposite limiting case, (3 — > 0, both quantities, r m and r 2 , are approaching zero. 
One can thus conclude that the rotationally symmetric self-trapped quantum probability 
density is converging toward a delta function, 

Bm/>(q ;) 0)-<5(q). (23) 

In Fig. HI we confirm the above conclusion numerically by plotting the radial profile of the 
quantum probability densities for decreasing small values of (3 — 10~ 4 , 10 -5 , 5 x 10~ 6 , 10 -6 . 
One can also see that for small (3 the quantum probability density is approximately taking 
a form of a rectangular function. 

Next, let us write the coupled dynamical equations of (jSJ) in polar coordinate, {r, 9}, to 
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FIG. 4: The profile of p(r) for various small values of /3. 

give us [3] 

m (—r- — ru 2 \ = —d r U, 
\dt J 

d t p+d r (pVr) =0. (24) 

where uo = d6/dt is the angular velocity and v r = dr/dt. In the last line we have used the 
fact that our dynamical system is rotationally symmetric. Let us now impose a stationary 
condition v r — 0. First, from the up most equation, the angular velocity is related to the 
quantum potential as 

io = y/d r U/(rm). (25) 

Hence, the angular velocity depends only on the distance, to = u(r). From the middle 
equation, since dgU = 0, one gets du/dt = 0, such that the angular velocity is constant. 
Finally, from the lower most equation, one gets dp/dt = dtp = 0. Hence, p(r; t) is stationary 
and spinning around its center with a stationary angular velocity field, uj{r). 

Going back to Cartesian coordinate, the rotating velocity vector field can be obtained as 

v = uj{t) x q = {— u(r)y, oj{r)x}. (26) 

This velocity vector field is divergence-less 
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This fact guarantees the existence of a twice differential) le scalar- valued function 99 (q), a 
stream function [17], such that 

v = {dytp, -d x tp}. (28) 

This means that the tangent of a curve <^(q) = ipo, where <po is arbitrary constant, is parallel 
to the velocity vector v. The set of curves y2(q) = tpo with various values of (po then make 
concentric circles with different radius all enclosing q = 0. See Fig. [5j 

If one further defines the quantum phase corresponding to this rotational flow as in Eq. 
dlj, one will have the following differential equations 

d x S = mdyip, d y S = —md x tp. (29) 

From Eqs. (T281) and fT2"9l . one gets 

dqS ■ <9 q y2 = m{v x (-v y ) + v y v x ) = 0. (30) 

Hence all curves S(q) = So are orthogonal to all curves y(q) = <fo, where So and ipo are 
arbitrary constants. Since the curves y?(q) = y^o are a set of concentric circles with different 
radius depending on the value of <po, then the set of curves S(q) = So are line rays emanating 
from the origin q = 0. See Fig. This means that at q = 0, S(q) is not uniquely defined. 
Hence, at q = 0, we have a phase defect or singular point such that Eq. is not valid. 

Finally, let us notice that Eq. (|29j) leads to an important fact that 5*(q) satisfies Laplace 
equation 

d 2 S = md x d y (f — md y d x (f = 0. (31) 

Hence, the quantum phase S(q) of our stationary wave function is a two dimensional har- 
monic function. On the other hand, the stream function satisfies the following Poisson 
equation 

d 2 qV = -9 q xv/0. (32) 

Keeping in mind all the above facts about the quantum phase, let us calculate the average 
quantum mechanical energy. Writing the wave function in polar form, ip = Rexp(iS/h), 
one obtains 

(E)= [ dqR( - —dl)R-i— [ dqRd Q R-d Q S 
Jm V 2m V mJ M 

+— / dq R 2 d Q S ■ d Q S - i— [ dqi)*d 2 S. (33) 
2m J M 2m J M 
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FIG. 5: The set of concentric circles are the curves of y(q) = tpo, whereas the set of lines emanating 
from q = are the curves of <S(q) = So- <^o and Sq are arbitrary constants. 



Notice that we have made clear that all the spatial integration is taken over a finite support 
M. inside which -R(q) 7^ 0. The first term on the right hand side is but equal to the average 
of quantum potential 

J M dqR (~ = Sm ^ P(q)f/(q) = °- (34) 

Next, recalling the fact that R depends only on r, and Eqs. (j4]), (1261) . one has 

dqR ■ dqS = —m—u)d r R + m—ud r R = 0. (35) 

Hence, the second term on the right hand side of Eq. (1331) is vanishing. Further, using again 
Eq. (j4j), the third term on the right hand side of Eq. (|33|) can be rewritten as 



1 

2m 



— I dq R 2 d q S ■ <9 q S 



J dq p^-mv • v 



M J M 

= [ dqp(ci)K{q)=K. (36) 
J M 

Here, we have defined a new quantity K = (l/2)mv 2 whose averaged over p(q) is denoted 
by K. Namely, K can be interpreted as the average kinetic energy of the Madelung fluid. 
Notice that since v(r) depends on U(r), so will K. We shall show later that K depends only 
on the value of /3, K = K(j3). Since (3 depends on U, one then conclude that K depends on 
U . Finally, by Eq. fl3TT) . the last term on the right hand side of Eq. (133]) is vanishing. In 
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total, we thus have the following decomposition of average quantum mechanical energy into 
the average quantum potential and average kinetic energy 



All the above facts tell us that the self-trapped, rotationally symmetric, spinning yet 
stationary quantum probability density we developed in this section comprises the class of 
wave functions that maximizes the Shannon information entropy, given the average quantum 
mechanical energy of the free particle. It is thus the most likely wave function of a single 
free particle in two dimensional space with quantum mechanical energy (E) = U + K. Of 
course we have to remind ourselves that we have worked in the frame work of Madelung 
fluid which we assume as the generalization of quantum mechanics allowing rotational flow. 

IV. A CLASS OF SELF-TRAPPED, NON-SPINNING YET STATIONARY WAVE 
FUNCTIONS 

Let us discuss the behavior of U and K of our stationary-spinning wave function as we 
vary {3. From the stationary condition of Eq. (1251) . it is clear that for finite value of /3, the 
rotating velocity of the Madelung fluid is infinite at the boundary of the support, dAi. Yet, 
one should again keep in mind that along this circle of radius r m , r 2 = x 2 + y 2 = r^, the 
quantum probability density is vanishing. For our stationary state, one can then show that 
the average kinetic energy is finite for non-vanishing value of /3 as follows 



where in the second equality we have used Eq. (I25I) . in the third equality we have used Eq. 
ffTU]) and in the fourth equality we used the partial integration. We have thus an important 
result which shows that the average kinetic energy K depends only on /3, inversely. For non- 
vanishing (3, assuming a system with finite average quantum mechanical energy (E), then 
the average quantum potential, U = (E) —m//3, is also finite. Moreover, given a finite value 
of average quantum mechanical energy, the definite positivity of U leads to the existence of 
a lower bound for the possible values of (3 as (3 > mj (E), such that r m is also non- vanishing. 



(E) = U + K. 



(37) 




(38) 
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FIG. 6: The average quantum potential, U, and average kinetic energy K versus j3. U is obtained 
for the solution of the differential equation (|12p with fixed boundary conditions, while K is obtained 
fromEq. (1381) . 

Fig. [6] shows the numerical plots of the average kinetic energy K and the average quantum 
potential U againstf (3 for the stationary-spinning wave functions which satisfy the differ- 
ential equation (fl~2l) with fixed boundary conditions and the stationary condition given by 
Eq. ( |38l) . One can see that, like K, U decreases as we increase the value of (3. Two limiting 
cases are of great interest. First is the case when /3 is vanishing, f3 — *■ 0. As discussed in the 
previous section, in this case the quantum probability density is approaching a delta function 
with vanishing support. From Eq. fl38l) one has lim^o-K" ~~ > °°- Numerical simulation also 
shows that the average quantum potential is approaching infinity, linig_>o U — > oo. Hence, 
this limiting case is not a good model or irrelevant for a single free particle of finite average 
quantum mechanical energy. 

Let us proceed to discuss the other extreme limiting case of infinite /3, /3 — > oo. Again, as 
discussed in the previous section, in this case, the radius of the support r m and the second 
moment r 2 of the quantum probability density are converging toward some certain finite 
values. This shows that the stationary quantum probability density and its corresponding 
quantum potential are converging toward some functions 

lim p(q;P) = Poo(q), 

lim U(q-f3) = U oc (q). (39) 

p— >oo 
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This also tells us that the average quantum potential, U, is converging toward a finite 
value, lim^oo U = Use, which is confirmed by numerical simulation. See Fig. [6j On the 
other hand, from Eq. (1381) . in this extreme case, the average kinetic energy is vanishing, 
lim^^oo K — > 0. Hence in this limiting case, the average quantum energy is exactly equal to 
the average quantum potential and thus finite 

lim (E) = lim U = U^. (40) 

f3— >oo /3^oo 

Moreover, since the average kinetic energy is vanishing, the quantum probability density is 
no more spinning. This situation is of great interest physically, since we have a lump of 
self-trapped quantum probability density which is non-spinning yet is still stationary. 

Let us investigate this last situation in more detail. Fig. [7J shows the profile of the 
quantum potential at increasing values of large /3, j3 = 1,5, 10, 100. One can see that the 
profile of quantum potential is approaching the form of a cylindrical tube with infinite wall, 
as we increase j3. Hence, for infinite value of /3, the quantum potential is finite and constant 
inside the support M., and infinite at dM.. Inside the support A4, the quantum force is 
thus vanishing. On the other hand, along the boundary dM, the quantum force is infinite. 
Yet one should keep in mind that along this circle boundary 8A4, the quantum probability 
density is vanishing. Hence, inside Ai, both the quantum force and the centrifugal force 
are vanishing. This ensures that Poo(q) is a non-spinning-stationary solution of the two 
dimensional Madelung fluid. 

Let us now proceed to discuss the profile of the wave function at the limit f3 — > oo. As 
discussed in the previous paragraph, in this limit the quantum potential inside the support 
is constant. Let us denotes the constant value of U inside M. by U c (See Fig. [7]). One 
therefore has 

(E) = £4o = J dq Poo (q)U c = U c . (41) 
Next, recalling the definition of quantum potential given in Eq. (j3J), inside A4, one has 

- |^*iUq) = t^oc(q) = (E)R 00 (q), (42) 

where we have denoted the quantum amplitude at infinite (3 by lim^oo i?(q; 0) = R^q). 
The above differential equation must be subjected to the boundary condition that along the 
boundary line of the support, dAi, the quantum amplitude is vanishing. Moreover, since in 
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FIG. 7: The radial profile of quantum potential for various large values of (3. 

the limit (3 — > oo the phase S is constant inside the support A4, multiplying both sides with 
exp(iS/h), one gets 

-|^^oo(q) = Woo(q), (43) 

where we have denoted ^oo(q) = -Roo(q) exp(iS/h). 

Eq. (|43p with the boundary condition described above is nothing but the stationary 
(time-independent) Schrodinger equation for a particle trapped inside a cylindrical tube 
classical potential whose bottom is flat and boundary is infinitely high. In polar coordinate, 
recalling the wave function is rotationally symmetric, one obtains the following differential 
equation 

^oo(r) + -d r ^(r) = -^P^co(r). (44) 
r nr 

The solution of which is given by 

^oo(r) = AS c (kr) exp(iS /h), (45) 

where S c (kr) = sin(Ax)/r is the Sine function, k = \j2m(E) /h 2 , A is a normalization con- 
stant and So is an arbitrary phase constant, k is thus a quantity of momentum dimensional. 
Fig. [8] shows the quantum probability densities for various increasing large values of j3 
obtained by solving Eq. (|T2|) with fixed boundary conditions, and its limiting case Poo( r ) 
given by Eq. (1431) . One can see that as (3 is getting larger, p(r; 0) is converging toward 
Poo(r). Already at f3 — 50, p(r; 0) is difficult to be distinguished from Pao(r). Next, from 
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FIG. 8: /j(q; (3) for various values of f3 which solve Eq. (|12h with fixed boundary conditions and 
its limiting case for infinite (5, Poo(r), given by Eq. (|45p . 

the boundary condition that the quantum probability density must be vanishing along the 
boundary line of the support, the Sine function S c (kr) must reach its first zero point at 
r = lim^oo r m (/3) = r^. One therefore has the relation between the average quantum 
mechanical energy of the single free particle and the radius of the wave function as kr^ = it. 

V. CONCLUSION AND DISCUSSION 

To conclude, we have first assumed the Madelung fluid dynamics as a generalization of 
Schrodinger equation for a single free particle moving in two dimensional space, allowing a 
rotational flow. We then ask for a class of wave functions whose quantum probability den- 
sity maximizes Shannon information entropy given the average value of quantum mechanical 
energy. We showed that there exist a class of solutions in which the quantum probability 
density is self-trapped, spatially localized with finite support, rotationally symmetric, spin- 
ning around its center and yet stationary. Hence it is a quantum vortex in xy— space. The 
stationarity comes from the balance between the attractive quantum force field of the trap- 
ping quantum potential generated by the quantum amplitude and the fictitious repulsive 
centrifugal force field of the rotating velocity field generated by the quantum phase. 

We then proceeded to show an asymptotic class of solutions in which the quantum proba- 
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bility density is no more spinning and yet still stationary. In this case, the quantum potential 
takes the form of a cylindrical tube with flat bottom and infinitely high wall. Moreover, 
we also showed that the wave function turns out to be the lowest eigenstate of the time- 
independent Schrodinger equation for a single particle trapped inside a cylindrical tube 
classical potential. In this sense, one can consider the class of spinning-stationary solutions 
of the Madelung fluid dynamics as the generalization of the notion of stationary states of 
the ordinary quantum mechanics. Yet, one should recall that in constrast to the ordinary 
quantum mechanics, we did not assume any classical potential to trap the particle. 

Some immediate interesting problems are then ready to be mentioned for future works. 
It is first interesting to study the case when the initial quantum phase does not satisfy 
the stationarity condition given in Eq. In particular, it is instructive to investigate 

the stability of the stationary solutions with respect to a small perturbation either on the 
quantum amplitude or quantum phase. In this context, it is interesting to compare the 
class of localized, stationary and spinning solutions developed in this paper with other 
localized wave phenomena like soliton etc. The case of three spatial dimension is then of 
great interest. Moreover, one is tempted to recover the quantization of energy in ordinary 
quantum mechanics using the formalism developed in this work. 
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